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SUMMARY 


In  many  equipments,  boiling  occurs  in  confined  space  such  as  the 
clearance  between  the  tube  and  the  support  plate  of  steam  generators. 
Corrosive  concentration  builds  up  at  the  boundary  of  dryout  zone  and  in¬ 
duces  severe  damage.  The  knowledge  on  this  kind  of  boiling  phenomena  is 
very  limited.  It  is  the  objective  of  this  research  to  understand  this 
fundamental  heat  transfer  of  this  problem  through  systematical  analysis 
and  experimental  studies. 

This  report  describes  mainly  the  analysis  performed  in  the  past 
year  at  Carnegie-Mellon  University.  The  forced  convective  two  phase 
flow  and  boiling  heat  transfer  In  confined  space  is  studied  using  the 
subchannel  analysis  which  is  an  approximated  method  but  effective  in  two 
phase  calculations.  Reasonable  prediction  of  the  dryout  phenomena  has 
been  achieved.  A  separate  analysis  is  performed  for  single  phase  flow 
heat  transfer.  This  differential  analysis  is  formulated  from  the  lubri¬ 
cation  theory  and  calculations  are  performed  for  both  the  line-contact 
and  two-point-contact  configurations.  Finally,  the  progress  in  experimen¬ 
tal  research  is  also  reported  here. 


CONCLUSION 

Both  the  subchannel  analysis  and  the  differential  analysis  calculate 
the  flow  field  in  confined  space  satisfactorily.  At  single  phase  flow 
the  differential  analysis  is  recommended.  At  two  phase  flow  the  sub¬ 
channel  analysis  is  very  effective. 

The  dry-out  pattern  and  the  temperature  field  in  the  confined  space 
have  been  predicted  by  the  subchannel  analysis  for  the  condition  of  line- 
contact.  The  results  are  reasonable  and  encouraging. 

The  fluid  flow  and  heat  transfer  at  single  phase  laminar  flow  have 
been  analyzed  by  the  differential  analysis.  The  formulation  and  the 
results  are  presented  in  non-dimensional  form  for  general  applications. 

In  the  coming  year,  the  experimental  data  will  be  obtained  to 
validate  the  analysis. 
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I.  INTRODUCTION 

Bolling  at  conventional  heat  transfer  surface  has  been  studied 
extensively  in  the  past  thirty  years  [1].  However,  at  the  same  time, 
the  boiling  in  confined  space  has  been  almost  completely  neglected. 

The  existing  information  [2]  [3],  is  limited  and  incomplete. 

In  many  equipments,  boiling  occurs  in  confined  space.  For  example, 
in  the  steam  generators,  clearance  is  allowed  between  the  heated  tube 
and  its  support-plate  where  the  tube  runs  through.  In  the  clearance, 
flow  is  reduced  and  boiling  occurs.  Another  example  occurs  in  nuclear 
reactor  core  where  the  structures  are  heated  by  nuclear  radiation. 

Boiling  occurs  at  the  narrow  space  among  the  structures. 

The  boiling  of  liquid  in  confined  space  leads  to  permanent  dryout. 
Corrosive  concentration  may  build  up  at  the  boundary  of  dryout  zone 
where  boiling  occurs.  Solid  deposits  may  also  occur  at  the  dryout 
boundary. 

The  understanding  of  the  thermal -hydraulic  related  corrosion  in 
confined  space  is  lacking.  This  is  primary  because  the  boiling  heat 
transfer  in  confined  space  is  not  really  known.  The  fundamental  under¬ 
standing  of  the  boiling  heat  transfer  in  confined  space  is  also  important 
to  the  advances  in  power  engineering  and  lubrication  engineering. 

In  view  of  this  need,  a  systematic  study  of  boiling  heat  transfer 
in  confined  space  is  performed  at  the  Dept,  of  Mechanical  Engineering 
of  Carnegi e-Mel  Ion  University  under  the  support  of  the  Office  of  Naval 
Research.  Both  experimental  and  analytical  research  are  performed. 
The  experimental  result  of  concentric  annulus  provides  the  fundamental 
information  on  boiling  heat  transfer  in  confined  space.  The  dryout 
pattern  in  eccentric  annulus  will  be  used  to  validate  the  prediction 
of  the  analysis. 


In  this  first  annual  report,  the  theoretical  analysis  and  its 
preliminary  results  will  be  presented  in  detail.  The  current  status 
of  the  test  loop  and  test  section  will  also  be  reported. 

II.  ANALYSIS 

Fluid  flow  and  heat  transfer  are  analyzed  for  the  narrow  con¬ 
fined  space  between  a  tube  and  the  tube  support  plate.  The  tube  can 
be  either  inclined  or  parallel  to  the  hole  of  the  support  plate.  When 
they  are  parallel,  the  annulus  may  be  concentric  or  eccentric.  Two 
practically  important  extreme  conditions  of  line-contact  and  two-point- 
contact  are  considered  as  shown  in  Figure  1.  Due  to  the  nature  of 
symmetry,  only  180°  of  the  total  confined  space  need  be  studied.  Add¬ 
itionally,  the  curvature  of  the  channel  can  be  neglected  in  analysis, 
such  that,  the  annulus  can  be  analyzed  as  a  flat  channel  with  varying 
channel  thickness.  The  typical  case  of  two-point-contact  is  shown  in 
Figure  2(a)  and  the  case  of  line-contact  in  Figure  2(b). 

Generally,  the  gap  thickness  h  for  an  eccentric  annulus,  see 
Figure  1(b),  can  be  expressed  acurately,  but  not  exactly,  by 

h  =  c  (1  -  e  cos  e)  (1) 

where  c  is  the  average  thickness,  e  is  e/c  where  the  e  is  defined 
in  Figure  1(b) . 


Therefore,  the  gap  thickness  h  is 

h  (xy)  *  c  ^1  -  cos  (f)j 

for  the  line-contact  condition  of  Figure  2(b), 

and  h  (xy)  *  c  £l  -  (2  ^  -  1)  cos(£)  j 


(2) 

(3) 


for  the  two-point-contact  condition  of  Figure  2(a). 
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A.  Subchannel  Analysis 

Viewing  the  flow  field  from  the  x-y  plan  as  many  individual 
subchannels  along  the  y  axis,  the  conservation  equations  can  be 
formulated  on  each  node  of  each  subchannel.  The  distinct  feature 
is  that  the  forms  of  the  conservation  equations  can  be  greatly 
simplified  as  compared  with  the  conventional  differential  analysis, 
such  that,  efficient  calculation  can  be  achieved  but  with  reasonable 
accuracy.  As  a  result  of  the  simplified  model,  the  subchannel 
analysis  is  able  to  calculate  the  complicated  two  phase  flow  and 
boiling  heat  transfer  effectively. 

In  the  present  analysis,  the  flow  field  is  calculated  on  the 
x-y  plan  with  the  varying  gap  thickness  h  considered.  The  velocity 
and  temperature  in  the  analysis  are  the  averaged  values  over  the 
thickness  h.  To  simplify  the  flow  field  calculation,  it  is  assumed 
that  the  axial  pressure  gradient  is  identical  at  any  point  at  a  same 
y  elevation  [4]. 


1 .  Formulation 
a.  Axial  Flow 

At  steady  state  two  phase  flow  in  the  confined  space, 

the  pressure  force  is  balanced  by  the  frictional  force.  That 

2 

means  -  Ap  =  *f,  \  f 

*  Dh  2  (4) 

where  the  f4  is  the  friction  factor  of  liquid  flow  with  a 

same  amount  of  mass  flow  rate,  G  is  the  mass  flux,  and  #  is  the 

two  phase  friction  multiplier.  Define  Fi  as  the  mass  flow 

rate  over  unit  length  in  the  x-y  plan  in  subchannel  i.  This  gives 
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Fi  a  p*“ihi 

Also,  the  pressure  gradient  at  the  elevation  j  is  the 
same  for  all  the  subchannels.  Set  that  as 

_  «  c2 

Ay  3 

From  the  above  three  equations,  the  local  axial  mass 
flow  becomes 


(5) 


(6) 


U 


(7) 


where  i  denotes  the  subchannel,  and  j  is  the  node 
notation  in  y  direction  as  illustrated  in  Figure  3. 

The  sum  of  the  axial  flows  should  equal  the  total 
mass  flow 


Assuming  the  two  phase  flow  in  the  narrow  gaps  as 
homogeneous,  the  frictional  multiplier  can  be  described  as 
$  =  1.0  if  x<0 

Cr)  ifx>0  (9) 

where  x  is  the  local  quality  to  be  evaluated  from  energy 
balance.  Better  forms  can  be  used  for  *  when  the  experi¬ 
mental  information  becomes  available  from  our  tests, 
b.  Cross  Flow 

The  cross  flow  between  subchannels  is  calculated 
from  mass  conservation.  Following  the  Figure  3(b),  the  mass 
conservation  gives 


FeAy  =  FwAy  +  Fsax  -  FnAx 


(10) 
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The  F$  and  Fn  are  known  from  the  axial  flow  calculations. 


The  boundary  condition  of  symmetry  at  x  ■  o  gives 


i x  *  o  w  (11) 

Therefore,  all  the  cross  flows  at  each  subchannel  can  be 
evaluated  easily  starting  from  the  locations  x  =  o  to 
x  *  iR. 


c.  Energy  Balance 

Following  the  illustration  of  Figure  3(b),  the 
energy  balance  for  a  control  volume  ij  can  be  written  as 


FsAxHs  +  F^AyHy,  +  q-j  jAxAy  =  FfiAy  H ^  +  F^xH^.  (12) 
where  H  is  the  enthalpy  and  q^-  is  the  heat  flux  from  walls. 
The  Fn  can  be  evaluated  from  equation  (10) 

FnAx  =  FwAy  +  FsAx  -  FeAy  (13) 

Substituting  the  Fn  into  equation  (12)  the  H-jj  can 
be  derived 


Hij  =  Hs  +  ^ij  ^/Fs 


+  F^Ay 
Fsax 


(14) 


+  FeAy 

The  third  term  at  right  hand  side  of  this  equation  will  be 
set  to  zero  if  Fw  is  flowing  in  the  negative  x  direction; 
the  fourth  term  will  be  set  to  zero  if  Fe  is  flowing  in 
positive  x  direction. 
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With  the  known  enthalpy  the  local  quality  can  be 
evaluated. 

*ij  =  (Hij  -  AHsub)  /  Hfg  for  0<X<1  (15) 

where  ^HsujJ  is  the  enthalpy  of  inlet  subcooling. 

2.  Numerical  Method 

The  calculational  procedure  is  shown  in  Figures 
4  and  5.  Iteration  is  performed  to  achieve  the  conver¬ 
gence  of  flow  rate  on  the  line  j,  the  pressure  drop 
across  the  confined  space,  and  the  local  quality  at  all 
the  nodes.  Due  to  the  simplicity  of  the  model  and  formu¬ 
lations,  the  convergence  is  rapid. 

3.  Results 

A  computer  program  has  been  developed  for  the  sub¬ 
channel  analysis  of  two  phase  boiling  in  confined  space. 
The  condition  of  line-contact  (refer  to  Figure  2(b))  is 
calculated  for  two  phase  boiling  with  dry-out. 
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The  two  phase  flow  and  boiling  are  calculated  for 
condition  of  line  contact  at  constant  heat  flux  condition. 

The  stream  line  of  the  fluid  flow  is  shown  in  Figure  6. 

Although  dryout  appears  at  locations  near  the  line  contact, 
since  the  gaps  there  is  small  and  the  dryout  pattern  is 
almost  parallel  to  y  axis,  the.  stream  lines  are  generally 
not  deviated  appreciably  from  a  straight  line. 

The  constant  temperature  contour  and  the  dryout  pattern  at 
the  line  contact  condition  is  shown  in  Figure  7.  At  the 
location  near  180°,  the  gap  is  large,  the  flow  is  strong, 
and  the  fluid  temperature  does  not  rise  much  from  the  inlet 
value.  Near  the  line  of  contact,  the  fluid  temperature  rises 
drastically  and  reaches  dryout  in  short  distance  of  flow. 

The  dryout  region  expands  rapidly  then  becomes  parallel  to 
the  y  axis.  It  is  believed  that  at  higher  system  pressures, 
the  boundary  between  liquid  region  and  dryout  region  will  be 
a  distinct  zone  of  two  phase  flow. 

Finally,  it  is  important  to  point  out  that  the  above 

analysis  contains  many  assumptions.  The  justfication  of 

these  assumptions,  or  the  future  modification  of  the  model, 

will  be  based  upon  the  comparison  with  experimental  data. 

However,  when  the  model  becomes  mature,  the  boiling  in  con¬ 
fined  space  can  be  predicted  for  any  fluid  at  any  working 

conditions. 


B.  Differential  Analysis 
1.  Formulation 

Precised  conservation  equations  can  be  established  in  the 
differential  form  for  the  fluid.  Due  to  the  nature  of  the 
narrow  flow  channel,  approximations  following  the  lubrication 
theory  will  be  adopted.  Although  the  present  analysis  is 
limited  to  the  single  phase  fluid  flow,  the  results  could  be 
used  as  a  reference  to  validate  the  more  approximated  sub¬ 
channel  analysis  which  has  been  used  by  us  to  study  the  two 
phase  boiling  phenomenon  in  the  confined  space. 

If  the  fluid  is  incompressible,  the  mass  conservation 
equation  becomes 


aw  .  Q 
ax  ay  3z 


or 


aw 
a  z 


3U_ 

ax 


3 v 

ay 


(16) 


(17) 


Themomentum  equations  can  be  greatly  simplified  according 
to  the  order  of  magnitude  analysis  considering  the  narrow 

gaps.  At  laminar  flow  they  become 

2 


ap 

3X 

a  u 

a  U  - 

*  3Z2 

(18) 

IE 

ay 

32V 

*  ^ 

(19) 

ap 

_  A 

3Z 

*  0 

(20) 

Parabolic  velocity  profiles  can  be  obtained  through  the 
integration  of  equations  (18)  and  (19).  Then,  the  velocity 
profiles  can  be  substituted  into  equation  (17)  and  integrate 


I 
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the  dz  from  0  to  h.  The  resulting  equation  becomes 

a  f"3  !£.■)  ,  _»/h»  o 

3X  V  3 x  /  3y\  3y/ 

which  is  the  Reynolds  equation  of  lubrication  theory. 


In  this  equation,  both  h  and  p  are  function  of  space. 
The  boundary  condition  for  forced  flow  is 


P  (x,o)  =  P0 

(22) 

p  (x,Ly)  =  o 

(23) 

—  (o,y)  =  o 

3X 

(24) 

3p 

—  (wR.y)  =  o 

3X 

(25) 

Once 

the  pressure  field  p(xy)  has  been  solved,  the 

"averaged* 

'  fluid  velocity  in  the  gap 

can  be  evaluated. 

u  =  - 

IZy  3X 

(26) 

v  =  - 

-1  IE  h2 

12u  3y  n 

(27) 

This 

set  of  equations  can  be  non 

-dimensional ized 

by  using 

x  =  “1 
ttR 

(28) 

Ly 

(29) 

h-  il 

c 

(30) 

p  =  £ 

Po 

(31) 

y  _  uit  itR 

PoC^ 

(32) 

V-  ^ 

(33) 

and 

0 

S  *  Ly 

(34) 

irR 

v* 
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We  can  further  define 

Pe  =  Re*PraVoLy/a  (35) 

where 

V0  =  PoC2/yLy 

Finally  the  governing  equation  for  pressure  field  becomes 


if  <*3  if  >  +  if  ^  17>  ’  0  (36) 

where  h  =  [1  +  (1  -2Y)]cos(irx)  for  2-point  contact  with 

(371 

boundary  conditions  ' 

P(X,0)  =  1  (38) 

P(X,1)  =  0  (39) 

|f  (0,Y)  *  0  (40) 

|f  ( 1  * Y )  -  0  (41) 


From  the  solved  pressure  field,  the  velocity  field  can  be 
found . 


U  =  "12  aX  (42^ 

V  =  -  -J2  if  (43^ 


The  energy  equation  can  be  written  in  terms  of  gap  averaged 
conditions.  The  final  energy  equation  becomes 


S2  U 


II 

ax 


v  II 
¥  3Y 


h  Pe 


(44) 


(This  equation  only  for  constant  heat  flux.) 
C  2 

where  Pe*  =  Pe*  (^p-) 


where 


T  = 


k(T-To) 

q  C 
w 


(45) 


-  > 
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with  initial  conditions 

T  =  0  at  Y  =  0  (46) 

|y  *  0  at  X  =  0  (47) 

2.  Numerical  Method 

The  above  formulation  have  been  written  into  finite 
difference  form  to  solve  for  the  pressure  field.  Different 
x  and  y  increments  are  allowed  in  the  program.  The  whole 
equations  are  solved  by  Alternative  Directional  Implicite 
method  [5]  to  speed  up  the  convergence  of  the  solution. 

Due  to  the  special  boundary  condition  at  the  line  x  =  o  and 
x  =  uR  (equations  (24)  and  (25)),  the  implicite  calculation 
is  not  performed  for  these  two  lines  during  the  sweep  of 
A.D.I.  method.  During  the  iteration,  the  implicite  equations 
on  a  same  line  are  arranged  in  the  form  of  tri diagonal  matrix. 
Then  Gauss  elimination  method  is  performed  for  the  solution 
of  the  pressure  field. 


3.  Results 

It  is  obvious  that  the  flow  field  of  single  phase  flow 
at  the  condition  of  line  contact  will  have  straight  stream 
line  along  the  y  direction.  For  a  same  pressure  difference 
across  the  tube  support  plate  the  overall  flow  rate  of  this 
differential  analysis  is  about  the  same  value  as  that  of  the 
subchannel  analysis.  This  confirms  the  compatibility  of 
these  two  approaches. 


The  typical  flow  field  of  the  two-point-contact  condi¬ 
tion  is  shown  in  the  Figure  8.  The  arrows  indicate  the  flow 
direction  and  magnitude.  The  curves  show  the  axial  component. 
They  possess  the  feature  of  symmetry  and  indicate  the  turning 
of  the  flow. 

The  temperature  field  of  two-point-contact  condition  at 
constant  heat  flux  condition  is  shown  in  Figure  9  for  the 
case  of  S  =  1.11.  The  contacting  points  show  high  fluid 
temperature.  Near  the  contact  point  of  (1.0),  hot  fluid 
extends  to  downstream,  but  diminishes  gradually.  Near  the 
point  (0.1)  the  hot  region  is  rather  localized.  Viewing  the 
overall  temperature  distribution,  it  is  clear  that  the  dry¬ 
out  phenomena  will  not  be  as  severe  as  that  of  the  line- 
contact  conditions. 

EXPERIMENT 

A.  Test  Loop 

In  the  past  year,  much  effort  has  been  devoted  to  the 
design  and  construction  of  the  test  loop.  At  the  present 
time,  the  loop  is  completed.  Test  runs  indicate  that  the 
design  requirements  are  fully  satisfied. 

The  schematic  of  the  loop  is  shown  in  Figure  10  with 


the  loop  charactersistics  summarized  in  the  Table  1. 


VACUUM 

(VENT) 


Figure  10 
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TABLE  1.  LOOP  CHARACTERISTICS 

Working  fluid  F  -  113 

System  Pressure  0.1  M.O  MPa 

Volume  flow  rate  0^68  LPm 

Preheater  Power  12  KW 

In  order  to  perform  boiling  experiments  at  constant 

pressures,  an  accumulator  is  attached  to  the  loop  to 

absorb  the  volume  expansion  of  the  two  phase  fluid  during 

the  tests.  Large  flow  by-pass  is  installed  parallel  to  the 

test  section  to  maintain  the  flow  stability  in  the  test 

section. 

The  instrumentation  system  has  also  been  established. 

The  instrumentation  rack  includes  the  digital  temperature 
displays,  turbine  flow  meter  readings,  and  the  strip  chart 
recorder  for  loop  temperatures.  The  Hesse  pressure  gauges 
are  mounted  on  the  loop. 

B.  Test  Section 

The  test  sections  have  been  designed  and  they  are  under 
fabrication.  The  schematic  of  pool  boiling  set  up  for 
studying  the  tube  sheet  crevices,  is  shown  in  Figure  11. 

The  test  section  is  heated  directly  with  an  internal  moveable 
thermocouple  for  measuring  the  inner  wall  temperature.  Through 
heat  conduction  calculation,  the  outer  wall  temperature  can  be 
calculated  and  then  the  local  heat  transfer  coefficient  is 
obtained.  The  pool  can  be  pressurized  to  about  0.6  MPa 
working  with  F  -  113  or  water. 
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The  schematic  of  forced  convective  test  section  is 
shown  in  Figures  12  and  13.  The  quartz  test  section  is 
changeable  (one  of  2.5  cm  length,  and  three  of  7.5  cm 
length  with  diffetent  I.D.'s). 

When  the  tube  is  heated  under  a  constant  heat  flux 
condi ton,  the  study  of  concentric  annulus  will  be  performed. 
The  test  sections  can  also  be  heated  under  constant  temp¬ 
erature  condition  using  a  heat  pipe  for  uniform  temperature 
control.  At  this  condition,  the  experiments  of  line-contact 
eccentric  annulus  will  be  performed  and  the  dryout  pattern 
will  be  studied.  The  dryout  information  will  be  compared 
with  the  prediction  of  our  subchannel  analysis. 


IV.  WORK  TO  BE  DONE 

The  heat  transfer  analysis  will  be  continued  for  the  conditions 
of  constant  wall  temperature  conditions.  Systematic  study  of  the 
flow  and  heat  transfer  processes  under  various  conditions  will  be 
studied  through  these  analysis  to  achieve  a  better  understanding  of 
the  boiling  phenomena  in  confined  space.  Experimental  research  will 
be  the  major  effort  in  the  coming  year  and  the  data  will  be  used  to 
validate  the  subchannel  analysis. 
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V.  NOMENCLATURE 


u 

IT 

u 

v 

7 


average  gap  thickness  of  the  confined  space 
the  pressure  gradient  at  location  j 
hydraulic  diameter 
eccentricity 
frictional  factor 

mass  flow  rate  over  unit  length  in  x-y  plan 
mass  flux 

local  gap  thickness  of  the  confined  space 
non-dimensionlyzed  gap  thickness,  h/c 
enthalpy 

latent  heat  of  vaporization 

enthalpy  of  inlet  subcooling 

total  length  of  the  confined  space 

total  mass  flow 

local  pressure 

non-dimensional  pressure 

the  pressure  drop  across  the  confined  space 

wail  heat  flux 

radius  of  the  heated  tube 

aspect  ratio  of  the  confined  space,  Ly/SR 

local  axial  velocity 

local  axial  velocity  average  over  gap  thickness 
non-dimensional  axial  velocity 
local  velocity  in  x  direction 
local  velocity  averaged  over  gap  thickness 


V 


non-dimensional  velocity  in  x  direction 

wr  velocity  in  z  direction 

x  coordinate;  or  local  two  phase  quality 

X  non-dimensional  x  coordinate 

y  axial  coordinate 

Y  non-dimensional  y  coordinate 

z  coordinate 

Greek  Symbols 

a  thermal  dlffusivity 

e  non-dimensional  eccentricity,  e/c 

e  angle 

y  viscosity  of  the  liquid 

p  fluid  density 

*  two  phase  frictional  multiplier 
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